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CP I Stochastic quantisation normally involves the introduction of a fic- 

titious extra time parameter, which is taken to infinity so that the 
system evolves to an equilibrium state. 

qh| In the case of a locally supersymmetric theory, an interesting new 

possibility arises due to the existence of a Nicolai map. In this case 

bJ!)' it turns out that no additional time parameter is required, as the 

existence of the Nicolai map ensures that the same job can be done by 
KJi \ the existing time parameter after Euclideanisation. This provides the 

j_j ' quantum theory with a natural probabilistic interpretation, without 

^ any reference to the concept of an inner product or a Hilbert space 

structure. 



1 Introduction 

In 1981 Parisi and Wu proposed a radical new way to obtain Green's functions 
in Euclidean quantum field theory [Tj. Known as stochastic quantisation, this 
approach enjoys certain conceptual and practical advantages over more tra- 
ditional ones. In this paper, we show how this approach arises naturally for 
locally supersymmetric models and provides them with a natural probabilis- 
tic interpretation of the wave function without reference to the concept of an 
inner product. 

Stochastic quantisation normally requires the introduction of an extra 
parameter t, referred to as the fictitious time. Evolution with respect to t 
is assumed to be stochastic, governed by a type of Langevin equation which 
gradually forces the system into thermal equilibrium. The static probabil- 
ity distribution arising in the t — * oo limit can then be identified with the 
Euclidean path integral measure and used to calculate the Euclidean Green 
functions. 

In the case of a supersymmetric theory, an interesting new possibility 
arises due to the existence of a Nicolai map. (This is a transformation which 
converts any supersymmetric theory into a non-interacting bosonic one [2], 
and which generally has the form of a stochastic differential equation in 
Euclidean theories.) In this case, there is no need to introduce a fictitious 
time parameter; the presence of the Nicolia map means that the same role 
can be performed quite satisfactorily by the physical time parameter from 
the original theory. 

Regarding the Nicolai map as a stochastic process leads to an interpre- 
tation of certain components of the wave function (in an appropriate repre- 
sentation) as probability densities whose integral is conserved in Euclidean 
time. This provides a simple stochastic interpretation for supersymmetric 
quantum theories, without the need to define an inner product on the space 
of states or appealing to the concept of a Hilbert space. This is particularly 
useful in the context of quantum cosmology, where it is difficult to identify 
the inner product [3]. 

In general, finding Nicolai maps for supersymmetric theories is not an 
easy task. In fact, it may be easier to find the associated Fokker-Planck 
equation directly, thus by-passing the difficulties of explicitly constructing 
Nicolai maps. 

This paper is organised as follows: A review of stochastic quantisation 
of a scalar field is given in §2J For a supersymmetric theory, the notion 



of Nicolai map is introduced in §31 In §H the canonical formulation of a 
one-dimensional supersymmetric non-linear a model is developed. Super- 
symmetry is always broken when boundaries are present. In ^we will see 
that there is a boundary correction to the action that restores invariance un- 
der a sub-algebra of supersymmetry transformations. Derivation of Nicolai 
map and Fokker-Planck equation in case of a rigid supersymmetry is done in 
§H1 In this section, it also has been discussed that the result can be obtained 
directly from canonical quantisation. Section ^ is devoted to quantisation 
of a locally supersymmetric model. 

2 Stochastic quantisation of scalar field the- 
ories 

The main idea in stochastic quantisation is to view Euclidean quantum field 
theory as the equilibrium limit of a statistical system coupled to a thermal 
reservoir. This system is assumed to evolve with respect to a fictitious time 
variable t and approaches the equilibrium limit as t — *> oo. The coupling to a 
heat reservoir is simulated by means of a stochastic noise field which causes 
the original Euclidean field to wander randomly on its manifold. In the 
equilibrium limit stochastic averages become identical to ordinary Euclidean 
vacuum expectation values. 

In this section we outline the stochastic, quantisation of scalar field theory 
using the approach of Parisi and WupJ. For a more complete discussion, the 
reader is referred to the review article by Damgaard and Huffel [4j. 

The basic idea is as follows: 

i) We imagine that each field (j){x) depends on an additional coordinate, 
the fictitious time t 

0(x)^0(x,t). (1) 

ii) We suppose that the evolution of 0(x, t) with respect to the fictitious 
time t is described by a stochastic differential equation that allows for 
relaxation to equilibrium. Specifically, one postulates that the evolu- 
tion is governed by the Langevin equation 

dt ~ 50(x)+^^'^' ^' ^^ 



where 5'[0] is the Euchdean action and r]{x, t) is a Gaussian white noise 
with correlations given by 



(r/(x,t)), = 

(r/(xi, ti)r/(x2, h))^ = 2,5"(xi - X2W1 - ^2) (3) 

iii) Given some initial condition at t = to and some realisation of r7(t), equation 
([2]) has a unique solution 0^(x, t). The correlation functions of 0^ are 
then obtained as Gaussian averages over all possible realisations of ?7(t); 

{U^,M) ■ ■ ■ 0.(x., h)), = ^^^^-|j,n,,,,.(,,) (4) 

As t -^ 0, equilibrium is reached, and the (equal time) correlation 
functions of 0^ tend to the corresponding quantum green functions 

lim {(j){xi, t) ■ ■ ■ (j){xk, t))r, = {(p{xi) ■ ■ ■ 0(xfc)) (5) 

An alternative to the Langevin approach is to study the Fokker-Planck formu- 
lation. In this approach, the stochastic averages are represented as functional 
integrals 

(0(xi, t) ■ ■ ■ <j){xk, t)), = J D<j)f{<l), t)<j){x,) ■ ■ ■ 0(xfc) (6) 

where the probability density functional f[(j){x),t] is a solution to the Fokker- 
Planck equation 

^-^=rd-x^(^^+^]f (7) 

dt J S<P{x, t) \ 5(t){x, t) S(P{x,t)J^' ^ ^ 

The physical correlation functions are then calculated using the equilibrium 
probability density functional 

/l^l=lim/|^,i] = r'W = -^. (8) 

If we define the operator 

= - /rf-a:— ^f-#- + --^), (9) 

' ^-"(X,t) V(50(x,t) 50(x,t)/' ^ ^ 
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then the Fokker-Planck equation (jTj) can be rewritten simply as 

df{x,t) 



dt 



-Of{x,t). (10) 



One can show that the resuhs obtained from stochastic quantisation agree 
with those obtained by more conventional methods [1]. Now for a system 
with a single degree of freedom x, and a potential V{x). The Fokker-Planck 
equation reads 

3 Supersymmetric models and Nicolai maps 

Nicolai has shown [2J that there is a map which transforms any supersym- 
metric theory to a free bosonic theory. The simplest example of Nicolai map 
is the Langevin equation, although in general for a typical supersymmetric 
action the Nicolai map will not have such a simple form [3J. 

The Nicolai's theorem essentially states the following: Take a supersym- 
metric theory and integrate out the fermion fields in the path integral. This 
contributes a non trivial determinant factor to the (bosonic) path integral. 
There now exists a transformation of the bosonic fields whose Jacobian de- 
terminant exactly cancels the determinant of the fermion integrations, and 
which simultaneously transforms the remaining bosonic part of the action to 
that of a non-interacting bosonic theory. 

The transformation (Nicolai map) will be invertible if one impose an ap- 
propriate number of boundary conditions on the bosonic variables in the the- 
ory. However, imposition of such conditions break the supersymmetry since 
the supersymmetric variation of the Lagrangian produces a total divergence 
which yields a boundary term when integrated. 

Nicolai's theorem will not be applicable unless the supersymmetry algebra 
has some graded sub-algebra whose bosonic generators preserve the boundary 
[5]. Note that in the framework of quantum theory, boundaries are impor- 
tant since one is generally calculating transition amplitudes between specified 
boundary data. Hence, boundary effects cannot be neglected. 

In general, only the Euclidean version of a supersymmetric theory will ad- 
mit a Nicolai map which can be interpreted as stochastic differential equation 
describing the evolution of the system in Euclidean time. This suggests that 
in quantum theory one can interpret the wave function (in an appropriate 



representation) as a probability density function whose integral is conserved 
in Euclidean time. In fact, we are more interested in the static state which is 
finally reached, since it represents the ground state of the Euclidean theory. 

4 N=l supersymmetry 

Here we present the canonical formulation of a one-dimensional supersym- 
metric non-linear a model describing a particle moving in a curved configu- 
ration space. This example will be used throughout the paper to illustrate 
an approach which can be applied to quite general supersymmetric models. 
Suppose the position q(t) of the particle at time t is described by n 
coordinates g*(t). The action for the locally supersymmetric Euclidean model 
is 5 = / Ldt with the Lagrangian 

+NV,,,i?ij^ - Wig.jq^ + NV,) - x^,(g^ - Ng'%) 
+Nxxi^^ij' - ^iVi?,,H?V'¥V'- (12) 

where gij{c{) is the metric of the configuration space in which the particle 
moves, V{q) is a potential function defined on this space, and g* = dq^dt. 
The covariant time derivatives of the fermion fields are defined as 

D^P' = ij' + r.k^^q^ (13) 

D^, = %- P,,lp^q' (14) 

where V^jk is the usual symmetric Christoffel connection on the configuration 
space, compatible with the metric gij [6]. The Riemann curvature tensor is 

R^jki = dkT'^ji — diV^jk + V^jiT^k — T^^jkT'^rni- (15) 

In the Euclidean formulation described above, the lapse function A^ is a 
real-valued function of time t. The standard formulation can be obtained 
by making A^ imaginary, or equivalently by taking A^ = iN with A^ real. 
Then the Euclidean Lagrangian L also becomes imaginary and so it is nat- 
ural to describe the theory in terms of the real function L = iL, which is 
the Lagrangian for the standard formulation. It follows that the momenta 



and Hamiltonian in the standard formulation are related to their Euclidean 
counterparts (defined below) by identities of the form Hgtandard = i^EucUdean 
and H standard = iH Euclidean- Note also that in the standard formulation the 
Grassman variables ip^ and ip must be related by complex conjugation to en- 
sure unitarity in the quantum theory. (This is unnecessary in the Euclidean 
version, where unitarity is not a requirement.) 

Having clarified its relationship with the standard formulation of the 
model, we henceforth consider only the Euclidean formulation described 
above. 

There are two different ways to define the (Euclidean) momenta conjugate 
to the variables g*. If we regard the Lagrangian 012p as a function of the 

variables q^^ifj^^ip^ and velocities q^,ip^,ip^, then differentiating with respect 

to g* (while holding ip^jip^jip^ and ipi fixed) gives 



dq'' 



^ = N-'g,,-^ - xA - XV^. + Thii^ji^'^- (16) 



However, if we regard L as a function of g*, ?/'*, ?/'j, g*, and the covariant veloci- 
ties D%1)^ and Dip^ and differentiate with respect to g* (now holding ^*, ip^, Dip^ 
and Z)^j fixed) we instead obtain a set of covariantly defined momenta 

dL — 

= p,- Vh^ilj.^K (17) 



9g* 
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These prove to be more useful for our purposes. 

From (fT2l) . we find that the momenta conjugate to X; X ^iid A^ are not 
independent quantities, but are subject to the primary constraints 

H^ ^ H^ ^ Hjv ~ 0. (18) 

Since H^, H^^ and Hn vanish at all times, their time derivatives must also 
vanish. After applying the equations of motion, this requirement gives rise 
to the secondary constraints [7] 

Q^Q^H^Q (19) 

where we have defined the quantitie^j 

Q = ii^\lii + Vi), Q = t{li,-V0 (20) 



^the ordering of the terms is immaterial here, but becomes more significant in the 
quantum theory. Different orderings may be used, and correspond to different choices of 
measure on the configuration space. 



and 

H = -^g'^{U,U, - V,Vj) + V,,^^^ - ]^,,ki^FiP'ij''^ljK (21) 

The constraint functions Q, Q and Ti are first class, and therefore generate 
of gauge symmetries. In fact, Q and Q are the generators of supersymmetry 
transformations, while Ti is the generator of reparametrisations. The factors 
of i in the definitions fl20|) and the initial minus sign in the definition fl2Tl) 
are included so that Q, Q and Ti are identical to the constraints appearing 
in the conventional formulation of the model. 

The total Hamiltonian is then found to be just a linear combination of 
these constraints, and has the form 

H = -N{n + ixQ + ixQ) (22) 

where N , Xi ^-^^d x can be regarded now as Lagrange multipliers. It vanishes 
when the constraints (fT9|) are imposed. 

There are also second-class constraints relating ■?/'* and t/'j to their conju- 
gate momenta. These constraints are interpreted as strong equalities and can 
be used to eliminate these particular momenta from the theory. Of course, 
the presence of second-class constraints requires us to employ Dirac brackets 
rather than the more familiar Poisson brackets. 

Following [8j, the elementary Dirac brackets are found to be simply 

{q\ g^} = 0, {g', H,} = 5*„ {H,, H,} = i)^i^'R\,, (23) 

{#,7A^-} = 0, {^^,V^^.} = 0, {^\V^^.} = 5^,. (24) 

(25) 
from which we obtain 

{g,7^} = 0, {Q,7^} = 0. (26) 

and 

{Q,Q} = 2n. (27) 

These equations are the hallmark of supersymmetry. 
Defining the fermion number F = ip^ip^, we also have 

{F,^.}=^. {F,^^} = -^j (28) 



while 

{F,Q} = Q {F,Q} = -Q, {F,H} = (29) 

showing that supersymmetry transformations do not mix fermion number, 
and that fermion number is unaffected by reparametrisations. Moreover, 
fermion number is a constant of the motion, since 

{F, H} ^ 0. (30) 

We conclude this section by remarking that in the case of rigid supersym- 
metry, x(^)) xit) and N(t) are specified functions rather than Lagrange mul- 
tipliers. Hence the supersymmetry and reparametrisation generators Q, Q 
and 7i are not constrained to vanish, and nor is the total Hamiltonian (!22l) . 

5 Boundary corrections to the action 

Supersymmetry is always broken when boundaries are present, since the su- 
persymmetric variation of the Lagrangian is a total divergence which yields 
a boundary term when integrated. In this section we will see that there is 
a boundary correction to the action that restores invariance under a sub- 
algebra of supersymmetry transformations [9] . 

The classical trajectories are invariant under the following infinitesimal 
supersymmetry transformations: 



Sq' = 


ei:' + eiP 


(31) 


6^' = 


-eg''{Iii-Vi)-r,ki^^5q^ 


(32) 


#, = 


-e(n, + Vi) + P,fe^/g^ 


(33) 


Sx = 


N-H + 2txx 


(34) 


Sx = 


N'^Z + 2exx 


(35) 


SN = 


-2(ex + ex) 


(36) 



where e, e are anticommuting Grassmann variables. The variation of the 
action is given by the boundary term 



5S=- 
2 



Ir- 



e^,(iV-^g^ + g'^V^ - x^') + 'ei^^N-'g^.q' -V,- x^J , • (37) 



t2 



While it is customary to disregard boundary terms, we cannot afford to do so 
here. The supersymmetry transformation is broken by the boundary terms 
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(157|) . preventing us from using Nicolai's theorem. In order to be able to 
construct a Nicolai map, we must find a way to restore tlie invariance of the 
action under a sub-algebras of the supersymmetry generators [5] . 

Exact invariance of the action under a sub-algebra of supersymmetry 
generators can be restored if the Lagrangian is augmented or diminished by 
the total derivative 

LB = j^{y + \ii,ij^). (38) 

Indeed, under the infinitesimal supersymmetry transformations fl3m36p the 
variation of the integral 



t2 



Lb dt 



ti 



1_ .1*2 



is found to be 
1 



61 =- 
2 



e^N-'q^ + g^W, - W) - eij\N-'g,,q^ -V,- xA 



Consequently, the supersymmetric variations of the modified actions 

rt2 



S^ = S±I= f\L±LB)dt 
Jti 



are 



and 



6S4 



6S- 



ei;^{N~'cf + g'^V,-W) 



li^%N-'g,^q' -V,- xi'^) 



(39) 

(40) 

(41) 
(42) 



If we specify that e(ti) = ^(^2) = 0, then 55'+ will vanish exactly. In other 
words, the modified action 5*+ is exactly invariant under the subalgebra of 
infinitesimal supersymmetry transformations obtained by imposing Dirichlet 
boundary conditions on e. We will refer to this as the left-handed subalgebra. 

Similarly, the modified action S^ is exactly invariant under the subal- 
gebra of infinitesimal supersymmetry transformations obtained by imposing 
Dirichlet boundary conditions on e. We will refer to this as the right-handed 
subalgebra. 

Because the modified actions differ from the original only by boundary 
terms, the classical equations of motion are unaffected. (In fact, by adding 
a boundary correction to the action, we have really just performed a canon- 
ical transformation.) However, the different versions of the action give rise 
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to different expressions for the momenta. If we use tlie modified action 
5*+ = J{L + Lb) (wliicli is invariant under tlie left-lianded subalgebra of 
supersymmetry generators), tlien tlie covariant momentum conjugate to g* is 



^+ ^ diL + Lb) 






(9g* 

and we can write tlie supersymmetry generators as 

Q = i^'Uf, Q = i{Ut - 2V,)^. (43) 

On the other hand, if we use the modified action S- = J{L — Lb) (which is 
invariant under the right-handed subalgebra of supersymmetry generators), 
then the covariant momentum conjugate to g* is 



^_ ^ d{L - Lb) 






and we can write the supersymmetry generators as 

Q = tiP\nr + 2V,i), Q = tU-^. (44) 

Quantisation of the theory then involves representing the canonical vari- 
ables as operators and interpreting t. However before doing this let us briefly 
consider the special case of rigid supersymmetry. 

6 Rigid Supersymmetry and Nicolai Maps 

The local symmetry described above reduces to rigid supersymmetry if one 
fixes the values of the (non-dynamical) Lagrange multipliers as (for example) 

iV = l, x = 0, x = 0. (45) 

In order that these values are preserved under supersymmetry transforma- 
tions, the transformation parameters must be constant: 

e = 0, e = 0. (46) 

Because A^, x, X have fixed values and are not allowed to vary, they no 
longer act as Lagrange multipliers enforcing first-class constraints. Hence, 
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when the supersymmetry is rigid, the quantities Q, Q and Ti defined in equa- 
tions (120112 ip are no longer required to vanish. As a consequence, the Hamil- 
tonian may also be non-vanishing. 

The two invariant forms of the action now reduce to 

(47) 
and 



(48) 
According to Nicolai's theorem [2], the invariance of each of these actions 
under a subalgebra of supersymmetry generators ensures that the theory can 
be transformed into a free bosonic theory by integrating out the fermions. 
The transformation will generally have the form of a first-order differential 
equation, and so will only be invertible if initial conditions are imposed on 
the bosonic variables g*(t) in the interacting theory. However, for Nicolai's 
theorem to work, any initial conditions must be invariant under the same 
supersymmetry subalgebra as the action [5]. For the action S+, appropriate 
initial conditions are 

q\h) = ql, ^P\t,) = (49) 

since these (like 5"+ itself) are invariant under the left-handed subalgebra. 
For S- one must use instead the initial conditions 

q\h)=ql Mti) = o. (50) 

Integrating out the fermions (subject to the appropriate initial condi- 
tions), one finds that the weight of a particular path q{t) in the ensemble of 
possible paths is given by p^ 



'P±[Qit)] = |[^,^]exp(-i^±[g,^,^f 

= J±[q]expl^-^J^y^g,0±V^^){q^±V'^)^ (51) 

J±[g] =exp|±i-£[<7^^V^,, -i/?]}rft (52) 
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and -R(q) denotes the curvature scalar of the configuration space at the point 
q(t). In fact the functional J±[q\ turns out to be the Jacobian for the trans- 
formation ^"(t) I— > g*(t) defined by the differential equation 

^ = Tg-'V, + e\ ■ e (53) 

where e*a is a vielbein field on the configuration space with the property that 

9ije\e\ = Sab (54) 

and the differentials Ag* = dq^ + g^^T^jkdt and e*a ■ ^""dt are defined so that 
they transform covariantly in the Ito calculus. It follows by a change of 
variable that the weight for a given history ^° (t) of the new variable is 



V[m] = exp [-1 £' \SatCedt^ (55) 



1 ft2 1 

and so C,"'{t) can be interpreted as a white noise process with auto-correlation 



{nt)e{t')) = ^s{t-t') (56) 

that drives the motion of the particle in configuration space via the Langevin 
equation (J53l) . The probability density function f±{t,q) for the particle's 
position q at time t will then evolve according to the associated Fokker- 
Planck equation 



dt dq^ 



9"vm^v.,u 



(57) 



where g denotes the determinant of the matrix of components Qij of the 
configuration space metric. 

In fact this result can be obtained directly from canonical quantisa- 
tion. The evolution of the quantum state vector \^{t)) is governed by the 
Schrodinger equation, which in this case can be written 

||vl/(t)) = \H\^{t)) = -^QQint)) - ^QQMt)) (58) 

on account of the operator identity QQ + QQ = —2hH that follows from 
the form of the Dirac brackets and our choices of values for the Lagrange 
multipliers. 
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As remarked above, Nicolai's theorem is applicable only if the action and 
initial conditions are invariant under one of the supersymmetry subalgebras. 
Requiring invariance under the left-handed subalgebra generated by Q means 
that the initial state |'?/'(ti)) must be annihilated by the operator Q\ and since 
this operator commutes with the Hamiltonian it follows that Q\ilj{t)) =0 for 
all t and hence 



i\m) - -^wi*w> - i. 



«(*)> = -;;i3WI*W> = —^i'^tm - '^v.,W IfC*)) (59) 



where we have used fj43|) to write Q and Q in terms of the covariant momenta 
n^ associated with the Q- invariant action 5"+. In fact these momenta are 
naturally represented by — /iVj, where Vj denotes the covariant derivative 
with respect to q\ 

In general the wave function will have a number of components, corre- 
sponding to solutions with different fermion numbers. Each such component 
is naturally represented as a p-form on the configuration space, where n— p is 
the fermion number [II1[I2]- The fermion annihilation and creation operators 
ip'^ and ip^ then act on a form uo according to 

ip'^u) = dq^ A u, ipjU = hijU (60) 

where ijU denotes the contraction of the form u with the vector field d/dqK 
With the momenta n~ represented in the manner described above, it follows 
that 

ip'Iituj = -Muj, Uf^'uj = h^6uj (61) 

where d is the exterior derivative and 6 is it adjoint, the coderivative. So if 
the state |\I/(t)) is represented by the form u, then (15^ implies 

^ = -ld{hSuj - 2ivuj) (62) 

where iy denotes the contraction with the vector field g^^Vjd/dq'^. The Q- 
invariant initial condition (l49l) implies '?/'*|\I'(ti)) = and hence F|\l/(t)) = 
for all t (since the fermion number F = ipiip'^ is conserved). Consequently, 
uj will be an n-form, and at each point in the configuration space must be 
proportional to the elementary n-form: 

10 = f+{q)dq^ Adq"^ A-.-Adq". (63) 
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It is then easily verified that fl62p agrees precisely with equation (1571) for 

/+(q)- 

Alternatively, if one starts with the Q-invariant action S^ and initial 
conditions fl50l) then the Schrodinger equation fl58l) reduces to 



|i*(o> - -^QQim) - i. 



*(*)> = -7^2 WI*W> = ^^i^jn.-'/'V^Cn- + 2V,) i*^) (64) 



and so if the state is represented by a differential form uj then 



duj 
'di 



— = -lS{hdu-2dVAuj). (65) 



(Note that in this representation it is the momenta IIj , rather than that are 
n^, that are represented by the operator —hVi] hence ip''Il~u = —Tidoj and 
Ti'^ip 1^ = ^ 5^-) In this case, the initial conditions ip |^(ti)) = imply that 
u; is a 0-form; say uj = /-(q) and so this equation reduces to fl57|) with the 
minus sign chosen. The evolution equation fl62l) is then seen to be equivalent 
to(l57Dfor/_(q). 

Which of the two representations is more appropriate depends on the form 
of the configuration space potential ^(q). For definiteness, we suppose that 
the potential grows as |q| increases such that lim|q|_,oo V^(q)/ln |q| > 2hn. 
In this case, the Fokker-Planck equation (157|) admits a static normalisable 
solution 

/4q) = Aexp(^-^^ (66) 

that represents the limiting t ^ oo probability density function of a Brownian 
particle driven by the Langevin equation 

^ = -9''Vi +e\-C- (67) 

7 Quantisation with Local Super symmetry 

Our goal in this section is to quantise the locally supersymmetric model, 
and to interpret the resulting theory. However we start by reviewing the 
method described above for rigid supersymmetry (albeit in slightly different 
notation). 
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If the model has only rigid supersyminetry, quantisation means finding 
solutions of the Schrodinger equation 

±\m}='jrHmt)} (68) 

with H = —N{H + ixQ + ixQ) and fixed values of x, X ^.nd A^. For the 
Euclidean theory we choose A^ = 1 and x = X = ? ^iiid so the general 
solution has the form 

|^(t)) = e-*^/^|^o) (69) 

with |\E'o) an arbitrary initial state. 

Assuming that V{q) grows as |q| increases, we choose the initial state 
|\E'o) to be invariant under left-handed supersymmetry; Ql^&o) = 0. Because 
[H, Q] = 0, it follows that 

Ql^(t)) = Vt > 0. (70) 

As discussed in the last section, in this case (jHSD can then be represented as 
a Fokker-Planck equation describing the evolution of a conserved probability 
distribution. 

The static probability distribution can be obtained by taking the t — *> cx) 
limit. (This limit must exist, due to the assumed form of the potential V^(q).) 
In the current notation, this is represented by the state 

l^oo) = lim |^(t)). 

The time-independence of this state imphes that it has zero energy; 

^|^oo)=0. 

(Note that in this respect the Euclidean theory considered here differs from 
the conventional theory, which admits a ground state with a positive energy 
and hence an oscillating phase. However, because the Euclidean evolution 
operator e~*^/^ is hermitian rather than anti-hermitian, only a zero-energy 
state can have a constant norm.) 

Thanks to fl7U|l . this limiting state is also Q-invariant: 

Ql^oc) = 0. (71) 
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Moreover, since = {<il ^\2hn\^ oo) = (^oolQQI^oo) + (^oolQQI^oo) = + 
II (5|^oo)|P, this state must also be Q-invariant: 

Ql^oo) = 0. (72) 

This is readily verified in the Fokker-Planck representation considered earlier; 
one has only to show that the function /+(q) given by fl66l) is annihilated by 
the operator Q defined in f H3l) . 

The state I^E'oo) is thus annihilated by all the first-class constraints of the 
locally supersymmetric theory and thus represents an acceptable quantum 
state for this theory also. 

Let us summarise the argument so far. Thanks to Nicolai's theorem the 
supersymmetry of the model ensures that there exists a representation of the 
quantum theory in which the Euclidean Schrodinger equation takes the form 
of a Fokker-Planck equation governing the evolution of a conserved probabil- 
ity density function. By taking the (Euclidean) long-time limit one obtains a 
time-independent state satisfying all the constraints of the locally supersym- 
metric theory and represented by a static probability density function. This 
procedure is readily adapted (at least in principle) to quite general theories 
with local supersymmetry, and is naturally regarded as a type of stochastic 
quantisation. 

Apart from giving a procedure for finding quantum state, this approach 
also provides a natural probabilistic interpretation of this state space without 
reference to any particular inner product. Indeed, by reviewing the deriva- 
tion of the probability density function fl66l) . the reader can confirm that the 
inner product has not been used at all. This is potentially of importance 
for theories such as supergravity, in which the choice of the inner product is 
itself problematic. By using the method outlined above, one can (in prin- 
ciple) obtain not only the quantum state that satisfying all the constraints 
of the theory but also the corresponding probability density function on the 
configuration space. 



8 Discussion and conclusion 

In this paper, stochastic quantisation of a locally supersymmetric model is 
performed without introducing an extra fictitious time variable as is normally 
done in stochastic quantisation. 
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In the context of supersymmetric models, Nicolai maps exist for certain 
components of the wave function. This suggests an interpretation of the wave 
function as a probabihty density whose integrals is conserved in Euclidean 
time. In the quantum theory, the bosonic components of the wave func- 
tion, satisfy a type of Fokker-Planck equation and thus can be interpreted 
stochastically. However, since our supersymmetric theory is local one can 
by-pass the difficulties of explicitly constructing Nicolai maps by finding the 
associated Fokker-Planck equation directly via canonical quantisation with 
a stochastic interpretation. 

In conventional quantisation of a classical theory one has to construct 
the Hilbert space and define an inner product in the Hilbert space. The 
probability then is the squared modulus of the wave function. However, in 
the stochastic quantisation of locally supersymmetric theory, the solution of 
the Fokker-Planck equation is a probability density function. Thus, there is 
a ready-made probabilistic interpretation of the quantum theory without any 
reference to the concept of inner product as one see in the usual Hilbert space 
formulation of quantum theory. This should be a big advantage in quantum 
gravity and cosmology, where the choice of satisfactory inner product has 
long been viewed as one of the most fundamental problems. 
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